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1 Introduction 

The method of differential renormahzation (DR) has appeared recently ||l|. It is a 
method of renormahzation in coordinate space that yields directly finite amplitudes 
and leaves the space-time dimension unchanged. Thus it is expected to be especially 
suited for supersymmetric theories, where the dimension of the space-time appears to 
be crucial for preserving the symmetry. DR has already been successfully applied in 
several contexts: the Wess-Zumino model ^, lower-dimensional ||3[ and non-abelian Q 
gauge theories, two-loop QED [||, chiral models Q, non-relativistic anyon models Q 
and curved space-time [0]. Other formal aspects of the method have been developed 
in 1^, |l^ and different versions of DR can be found in |[TI[] . It is the purpose here 
to push even further the method and tackle a relatively complex problem such as the 
computation of {g — 2)i in supergravity (SUGRA). This is the first time that DR is 
applied to the calculation of a physical observable and to the calculation of gravitational 
corrections. As we will show our results preserve supersymmetry (SUSY) and abelian 
gauge invar iance. 

Although SUSY is not an exact symmetry of Nature, it is believed that any funda- 
mental theory must be originally supersymmetric. It is also known that when SUSY is 
made local it naturally includes gravity ||T^ . The resulting theory is non-renormalizable 
but it is constrained by the symmetries. One of the few finite calculations in gravity is 
the one-loop correction to the anomalous magnetic moment of the lepton {g — 2)i |l3| . 
In SUGRA not only is {g — 2)/ finite but SUSY requires it to be zero [0]. Therefore, 
the anomalous magnetic moment of the lepton besides being an observable, is also an 
ideal arena to check theoretical implications and to perform consistency tests of the 
methods of regularization in SUGRA. 

In a supersymmetric theory (g — 2)i vanishes because no such term appears in the 



Lagrangian of a chiral supermultiplet |14]. (This has been generalized to a set of sum 
rules valid for any charged supersymmetric multiplet [^, fB].) Hence, as long as SUSY 
is preserved, all quantum corrections must cancel order by order. Ferrara and Remiddi 
also proved explicitly that in global SUSY the one-loop QED corrections, order e^, do 
cancel. In this case two diagrams contribute (see Fig. 1), one exchanging a photon 
and the lepton, and the other a photino and the corresponding sleptons, and both are 
finite. 

The one-loop gravitational corrections are of order en"^ = SvreGiv, resulting from 
a graviton or gravitino exchange. Using dimensional regularization |jl^, Berends and 
Gastmans calculated the five diagrams where a graviton is exchanged (see Fig. 2) [|13|| . 
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Figure 1: Diagrams of order e^ contributing to {g — 2)i in superQED. 
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Figure 2: Diagrams of order ck? contributing to {g — 2); in SUGRA. A graviton is 
exchanged in diagrams D1-D5 and a gravitino in D6-D10. 



All five diagrams are infinite but their sum is finite. The finiteness of (^f — 2)i in a 
non-renormalizable theory such as gravity seemed miraculous. Del Aguila et al. [18| 
and Bellucci et al. ||l9| checked that when gravitation is embedded in a supersymmetric 
theory (unbroken), the contributions from the graviton and the gravitino cancel, as 
required by SUSY. Bellucci et al. also traced back to an effective chiral symmetry in 
the gravitino sector the finiteness of the gravitino contribution and then of the graviton 
contribution, if their sum has to vanish. Dimensional regularization does not yield a 
vanishing value for (g — 2)i. This is not so surprising for dimensional regularization 
is known to break SUSY. A (one-loop) SUSY preserving method such as dimensional 
reduction [^] is then required in order to obtain such a cancellation and this was shown 
to be the case. However, whether the finite contributions to {g — 2)i of the graviton 
and gravitino sectors are well-defined quantities for some unknown reason (symmetry) 



is a question which remained unanswered. The lack of manageable, alternative regu- 
larization methods preserving SUSY has prevented to address this question although 



the answer is expected to be negative as hinted in |l^, 21, 13]. 

In this paper we use DR to calculate the one-loop corrections to (g—2)i in unbroken 
SUGRA. We find that in DR the individual contributions of the graviton and gravitino 
sectors, although finite and opposite, are different from previous results. Thus, although 
SUSY is preserved, these two contributions are not well-defined separately, only their 
sum is. 

The method of DR is based on the observation that in coordinate space singularities 
arise when points in the quantum fields coincide. The idea is to write singular expres- 
sions as derivatives of less singular functions. One solves some differential equations and 
considers the solutions as a definition of the amplitudes in the sense of distributions, 
i.e., the derivatives are understood to act on test functions. In fact the method gives a 
definition of the expressions at the ill-defined points. The constants that naturally ap- 
pear in the differential equations play the role of renormalization scales. The resulting 
finite amplitudes satisfy the renormalization group equations. While retaining the spirit 
of the method, we have developed some aspects that may be considered unsatisfactory 
in the original version. For instance, gauge invariance had to be checked at each step of 
the computation and the renormalization scales chosen accordingly. In this paper we 
show that it is possible to preserve gauge invariance (no term proportional to the pho- 
ton momentum is generated at one loop) and supersymmetry {{g — 2)i has no one-loop 
corrections) if related singularities are treated always in the same way throughout the 
computations. Only one scale for each type of singularity is then needed. In particular, 
for the computation of {g — 2)i only the scale corresponding to logarithmic singulari- 
ties appears. Obviously it cancels out when all diagrams are summed up. These new 
features of the DR approach will be more extensively discussed in a simpler context in 
a forthcoming publication [ |22| . 

The plan of the paper is as follows. In Section Q the lagrangian describing a 
minimal superQED-SUGRA theory is presented. Section ^ is devoted to describing 
the constrained differential renormalization procedure which we will follow. Section B 
contains the detailed calculation of one diagram, where the main techniques used in 
the paper are applied, whereas in Section ^ we discuss briefly the calculation of the 
remaining diagrams contributing to (g — 2)i and give the final results. We finish with 
the conclusions in Section ^. Three Appendices gather the Feynman rules in Euclidean 
space and other technical details. 



2 The SUGRA Lagrangian 

The coupling of matter to supergravity has been extensively studied in the litera- 
ture [23|. We are interested here only in first-order gravitational corrections. We shall 
use the lagrangian of superQED-SUGRA obtained by imposing canonical kinetic terms 
(minimal Kahler potential and / function) and expanding the curved metric around 
the Minkowski one {rjap)- The interaction lagrangian in Minkowski space reads 

- eV2{\(t)\VL^ + h.c.) + (L ^ R)] , (2.1) 

Ceeg+eeg^ = -^/l"^[(i^(7a5/3 + T/?^^)^ + h.C.) 

-2e§(7,^^ + 7/3vl„)^] , (2.2) 

Ceeg+eeg-r = " ^[x''^l(^ ^ " "^)(/>J^7^,^ 

+ ex^PL 4^7^^ + h.c] + {L^R) , (2.3) 

+ i^-XY[^,4]x, + h.c.)], (2.4) 

^-» ^- <— 

where d= d— d, with d acting on the function on the left, and Vr^l = 2(-'^ ^ 75) 

are the chiral projectors. The kinetic terms are the canonical ones, with the photon 

in the Feynman gauge and the graviton in the de Bonder gauge. We use the following 

notation for particles and their corresponding fields: 



lepton 


-^ e , 


*, 


photon 


-^ 7 , 


-^fij 


graviton 


-^ 9 , 


t^lMU 



sleptons 


-^ eL,R , 


<pL,R 


photino 


-^7 , 


A, 


gravitino 


^9 , 


Xm- 



The lagrangian is written in Minkowski space for easier comparison with previous work. 
However we shall work in Euclidean space as is usual in DR. The corresponding Feyn- 
man rules are collected in Appendix A. 

3 Constrained differential renormalization 

The e^ and ck? corrections to the lepton-lepton-photon vertex in SUGRA are given 
by the diagrams in Figs. |^ and |^, respectively. Their expressions contain many differ- 
ent singular pieces. To obtain renormalized expressions, each singular piece must be 
substituted by a regular one according to the DR prescription. It is crucial, however, 
to perform all such renormalizations in a consistent way if the symmetries are to be 



preserved. In the literature this is taken care of by imposing certain relations among 
the renormalization scales introduced in every DR replacement. For instance, Ward 
identities in QED require that the 'logarithmic' and 'quadratic' renormalization scales 
in the one-loop vacuum polarization be equal, while the scales in the fermion selfenergy, 
Ms, and the vertex correction. My-, must satisfy the equality log j^ = \ iQ. 

Our approach here is to proceed in such a way that symmetries are automatically 
preserved. This can be accomplished if all singular pieces are written in terms of 
a minimal set of singular functions, which are renormalized afterwards. This ensures 
that related singularities are treated in an identical way, no matter in which diagram or 
position they appear. Within this scheme, differential renormalization should not break 
the symmetries of the bare theory. In Ref. [^2| it is shown, for instance, that with this 
approach the vertex Ward identity in spinorial QED and scalar QED is automatically 
satisfied to the one-loop level. Here, the final result will be directly compatible with 
both supersymmetry and U{\) gauge invariance. 

This 'constrained' differential renormalization relates diagrams of different topol- 
ogy. This is done 'separating' the points of diagrams with a smaller number of propa- 
gators. Let us classify the diagrams in Figs. IT] and in two classes according to their 
topology: 

i) Triangular diagrams: they are products of three propagators with at most four 
derivatives. Translation invariance and a systematic use of the Leibnitz rule allow 
to express them in terms of a set of functions defined as 

T[0]E^A(x)A(y)[0-'A(x-y)], (3.1) 

where x = xi — X3, y = X2 — x^ and Xi are the coordinates of the triangle 
vertices. A(x) is a scalar propagator and O^ is a differential operator acting on 
X of order < 4. Each of these propagators can be either massless or massive. 
(In the following, with the exception of Eq. (^^) below which is analogous to 
Eq. ( p.lp but for bubble diagrams, we reserve the notation A(x) for massless 
propagators.) In our calculation we find T-functions with one massless and two 
massive propagators, 

Ti[0] = A„(x)A„(y)[0"A(x - y)] , (3.2) 

and T-functions with two massless and one massive propagators, 

T2[0] ^ A(x)A(y)[O^A^(x - y)] , (3.3) 



where A(x) = 4^^ and /S.m{x) = _1^ ™- U™-^) fulfill the propagator equations 

n^A(x) = -5{x) , (3.4) 

(n--m2)A^(x) = -5{x) , (3.5) 

and Ki are modified Bessel functions |24]. The mass structure does not affect the 
singular behaviour. Hence, as far as renormalization is concerned, Ti and T2 can 
be considered identical (although some care is needed when both massless and 
massive DR identities are used, as discussed in Section |5|). 

ii) Bubble diagrams: they are products of two propagators and a delta function, 
and contain at most two derivatives. They can be written in terms of a set of func- 
tions whose general expression, omitting the index for any massive propagator, 
is 

By[0] = 5(y)A(x)0"A(x) , (3.6) 

and analogously for Bx, with x ^^ y. O^ acts on x and is now of order < 2. In our 
calculation B-functions always contain one massive and one massless propagators. 
Therefore the notation 

By[0] = 5{y)^{x)0^/\^{x) , (3.7) 

B,.[0] = ,5(x)A(y)OJ'A„(y) , (3.8) 

can be used in what follows. 

Although the structure of B-functions and T-functions is apparently different — and a 
priori they could be renormalized differently — , they are in fact related. B-functions 
can be expressed in terms of T-functions using the propagator equality ( |3.5| ). This 
substitution 'separates' points and allows to express a bubble diagram as a triangular 
one. In this way one is able to express all diagrams (and then all singularities) in terms 
of T-functions only. The procedure is described in detail in Section ^. There it is 
also proven that all the relevant singular T-functions can be written using one single 
singular function: T[n]. Therefore just one renormalization will be eventually required 
and only one arbitrary scale parameter will have to be introduced. Although the full 
correction to {g — 2)i is finite and does not depend on the renormalization scale, we 
calculate the contribution of each diagram separately in order to compare with previous 
results. 



4 A detailed example 

In this Section we present in detail the evaluation of the contribution to {g — 2)i of 

diagram D6 in Fig. ^. Using the Feynnian rules in Appendix A, this vertex correction 

reads (all indices, including the 7-niatrix ones, are in Euclidean space) 

2 
V^^\xi,X2,X3) = 2^[7„(^^i -m)PfiA^(xi -X3) 9^' A„(x3 - xa)] 

X [7^ r'""'7aA(x2 - xi)Pi(/' +mhfs] . (4.1) 

The factor 2 comes from the fact that the scalar particle propagating in the diagram 
can be either cl or cr. Translation invariance allows to write V^ as a function of 
X = xi — X3 and y = X2 — xs only: 

Vl,'\xuX2,X3) = i;l^)(x,y) 

= !^[d^^A{x-y)] 

X bai^'' - m.)PiJ7/37r7a 

x{^y + m)7/3(a;)(A^(x)A„(y))], (4.2) 

where dp = d^ — dp. The terms containing a 75 do not contribute to {g — 2)i and 
hereafter will be ignored. As explained in Section ^ all vertex corrections can be written 
in terms of derivatives of T-functions. Using the Leibnitz rule and the properties of 
the 7-matrices, we find 

ien^i^m" - 29^ • 9^)5; 7/3 - 2md~d-p)Ti[dp\ 

+ ((2m2 - Ad^ ■ dy)^p6^p - Amd-p6o^^ + 29; 9^ 7^ - 2md-6^p)Ty[do.dp\ 
+ (49^7q(5^t, - Am5a|35^^l + 2d~-^i35a'^)Ti[dadi3d^] 
+ 47^Ti[a^a^n]} . (4.3) 

The external derivatives must be understood in the sense of distribution theory accord- 
ing to the DR prescription: they act 'on the left' over test (wave) functions. Of the four 
T~functions in Eq. ([4. 3D, only Ti[9/3] is regular ( i.e., a tempered distribution). The 
other three are singular at 2; = y = and must be renormalized (regularized). This 
renormalization is done in three steps. 

First, the singular T-functions are split according to their tensor structure into 
trace and traceless parts: 

Ti[d^dp\ = Ti[d^dp--^5^pn] + -^5^pTi[u]^ (4.4) 



+ \{5o,pTi[d^u] + 5„^Ti[a^n] + 5p^Ti[do^u]) , (4.5) 

Tildo^dpU] = Ti[{dadp-^6apa)n] + ^6apTi[an]. (4.6) 

The traceless parts are two orders less singular. Thus, Tildadp—jdap^] and Ti[dadp d^- 
^{6apd^ + Sa-ydfs + Spyda)^] are finite and TiKdadp — i(^Q,^n)n] is 'logarithmically' 
singular. Ti[n] and Ti[nn] remain 'logarithmically' and 'quadratically' singular, re- 
spectively. Ti[5Q,n], which seems to be 'linearly' singular, is in fact 'logarithmically' 
singular because it can be written as a function of Ti[n]. Indeed, using the symmetry 
of T-functions under x <-> —y interchange (see Appendix B), 

Tild^D] = -^d-T,[a] . (4.7) 

Second, all singular T-functions are expressed as functions of Ti[n] and Ti[nn] 
only. Indeed, the other singular T-function left, Ti[{dadp — jdapO)^], can be written 



[-{d^dl + dldl) - -idldl + dld^p) + -5^p{d^ -dy-U-- D^)]Ti[D] 



+\j^Wa + di){di + di) - \5^p{d; + dy){d; + dy)] 



as is proven in Appendix B. Only the first term in the r.h.s. of Eq. ( [4.8D is singular. 

Third, Ti[n] and Ti[nn] are renormalized. The singularity of Ti[n] goes as the 
inverse of the distance to the origin to the fourth power and is easily renormalized. 



Inserting the propagator equation (3.4), one obtains 



Ti[n] = -[A^{x)A^{y)]6{x-y) 
1 m?Kf(mx 



(47r2)2 X 



2 



<5(x - y) , (4.9) 



and following Ref. |10| the renormalized T- function is 



rr^Rr 1 If/ Nrl/ . n.mKo(mx)Ki(mx) o, M^ ^, ,, ,,,^. 
Tf[°] = -7T^T2'^(^-y)[9(°-4"^ ) -^ '- + 7T^log—6{x)] , (4.10) 

where M = 2M/je is an arbitrary scale. The renormalized Ti[nn] can be found in 
Ref. |2^]. However it does not contribute to {g — 2)i because the corresponding term in 
Eq. ( [4.3| ) is proportional to 7^. Hence, the renormalization of this correction to {g — 2)i 
reduces to renormalizing Ti[n] as above. 

9 



Putting everything together we can now evaluate the contribution of D6 to the 
anomalous magnetic moment of the lepton. We substitute Eqs. {AA - ^^ ) and ( 4.1C| ) 



into Eq. ([4.3| ) ignoring the Ti[nn] term. Then we Fourier transform the renormalized 
vertex correction and put the lepton and the photon on their mass-shells. The part 
proportional to p^ — p'n, where p and p' are the incoming fermion momenta, gives the 
{g — 2)i contribution of this diagram. The Fourier transforms of the regular terms 
entering in (g — 2)/ are given in Appendix C. They add to 

(5-2)r = ^(-log4 + -)- (4.11) 

In other diagrams the proliferation of terms makes these manipulations rather cumber- 
some. However, a symbolic program developed for this purpose greatly simplifies the 
calculation. 

5 The complete calculation 

In this Section we describe the main aspects of the computation of {g — 2)/ for all the 
diagrams in Figs. || and |2[ 

Diagrams DA, DB, Dl and D6. The first two diagrams (Fig. |l|) are of order e^, 
and the rest of order bk^. In the example above we obtained the contribution of D6 to 
{g—2)i and showed how to treat expressions containing Ti-functions. In diagram Dl the 
same set of functions appears, so the procedure is completely analogous. In particular, 
quadratically singular terms are again proportional to 7^ and do not contribute to 
{g — 2)1. In fact, the same occurs for all diagrams because the quadratically singular 
functions are scalars containing the maximum number of derivatives (four in T's and 
two in B's), so the only Lorentz vectors left are Dirac gammas. Diagrams DA and DB 
are even simpler, for no singular terms contribute to (g — 2)i once the trace-traceless 
splitting has been performed. 

Diagrams D4, D5, D9 and DIO. These diagrams contain two massless and one 
massive propagators, and are written in terms of T2-functions. Here three derivatives 
appear at most, so there are no quadratically singular terms (not even in the 7^ part). 
All T2-functions are reduced to T2 [□] plus regular terms in a similar way as we reduced 
Ti-functions. Then the equation (|3.5| ) is used to write this singular function as 



T,p]=m'T,[l]-^^J{x-y), (5.1) 



which together with the DR identity |l| 



in 



allow to define the renormalized function: 

T^[a] = m'T,[l] + 4(4;^°^^^5(^ - y) • (5.3) 

This renormalization is the same as the one in Eq. (|4.10| ), so it involves the same mass 

scale M. Indeed, if the function in Eq. ( [4. 91 ) is expanded in the mass parameter, we 

obtain 

m?Kl{mx) 1 



r2 



+ n(m,x) , (5.4) 

x^ X"* 

where 7^(m,x) is of order m"^ and regular, and the r.h.s. of Eq. ( |5.4D can be also 



renormalized with the identity ( |5.2| ). In Eq. (4.10) we used a massive DR identity 
instead, in order to obtain more compact expressions, but we were careful to define 
the renormalized expression in such a way that it agrees with the one obtained from 
massless renormalization of the expansion in Eq. (5.4). Hence the scales in Eqs. ( 4.10| ) 
and (|5.3|) are the same. 



Diagrams D2, D3, D6 and D7. They can be written in terms of the functions 
By (D2, D6) and Bj. ( D3, D7). As discussed in Section y, one must express B- 
functions in terms of T-functions before renormalizing. This can be achieved using 
the equation ( p. 5] ) to 'separate' the delta function into a propagator, thus obtaining a 
triangular structure: 

By[l] = Siy) A{x)A^ix) 

= 6{y)A{x - y)Am{x) 

= -□^A„(y)A(x - y)A„(x) + m^Amiy)Aix - y)A^(x) 

= _(□?/ _ m2)Ti[l] - 29^Ti[5,] - Tip] . (5.5) 

Analogously we obtain: 

B,[5^] = -d-^{ay-m^)Ti[l] + {ny-m^)Ti[d,]-2dpyTi[d^] 

- d^Tip] + 2dyTi[d^d„] + Ti[5^n] , (5.6) 



Byid^d,] = -d'-d:{ay-m^)Ti[i] + d-;:{ay-m^)Ti[d,] + dUoy-m^)Ti[d^ 



- 2d^^d-,dyTi[d^] - 9^9^Ti[n] - {ay - m')Ti[d,d,] 

+ 2a^9^Ti[9,a,] + 2d^dyTi[d^d^] + d^Ti[d,a] + d^Ti[d^n] 

- 2dyTi[d^d,d^] - Ti[d^d,n] , (5.7) 

and similar formulae for Bx- Had we naively renormalized the B-functions indepen- 
dently, we would have a priori no control over the relation between the renormalization 

11 



scales of triangular diagrams Mt and of bubble diagrams Mb , and the choice Mt = M^ 
would break SUSY. 

Once all the relevant terms in a diagram have been renormalized, one only has to 
extract from the regular expressions the contribution to {g — 2)i. This is conveniently 
done by performing a Fourier transform and taking the appropriate on-shell limits, as 
described in Appendix C. 

Let us summarize the steps we have followed to evaluate {g — 2)^: 

1. All diagrams are written in terms of T- and B-functions. 

2. B-functions are rewritten in terms of T-functions. 

3. The singular parts of the T-functions are identified and renormalized. 

4. Finally, the contribution to {g — 2)i is extracted Fourier transforming, putting the 
external leptons on-shell and taking the g^ ^ limit. 

The explicit calculations have been carried out by hand and checked with a dedicated 
symbolic program. The necessary equalities are collected in the Appendices. 
The resulting contributions of the superQED diagrams are 

Note that they are finite, scale independent, and, of course, agree with previous results. 

The graviton and gravitino corrections to (g — 2)i are given in Table 1, together 
with the dimensional reduction and dimensional regularization results. In next Section 
we comment on the different contributions and compare them. 

At this point, as the total graviton and gravitino contributions are finite, one may 
wonder if it is possible to add all the contributions of the different diagrams to start 
with and work with non-singular expressions. Actually this requires treating all the 
diagrams in appropriate manner. Writing all the contributions in terms of T-functions 
our procedure allows to obtain a compact result free of singularities. We find for the 
part of the vertex amplitude contributing to (g — 2)i 

.2 

I [-2m^d~-fi3 + 4m35^^ - 8md^dp\ Ti[dp] 

-8md;,Ti[a] + 8m^d;,T2[i] 

+8md;,d^T2[dp] + 8md-T2p]} . (5.10) 

12 



-ZCK 



This expression is regular because the only singular terms are those proportional to 
Ti [□] and T2 [□], which have the same singular parts: 

T2[n] = Tip]-{A{x)-A„,{x))A„,{x)6{x-y) 

-(A(x) - Arnix))A{x)6ix -y) + m^Tsil], (5.11) 

and appear with opposite sign. 

Another important remark, further discussed in Appendix C, is the infrared diver- 
gent behaviour of the Fourier transforms of the T2-functions in the static limit q'^ -^ 0. 
It can be shown, however, that these divergencies cancel out in the sum Eq. (5.10). In 



fact we have checked that all diagrams, with the external leptons on-shell, are well- 
behaved in this limit. 

We have also verified the absence of gauge non-invariant terms, proportional to 
the photon momentum q^, when the leptons are on their mass shell. This corresponds 
in coordinate space to a symmetry under the x <-> —y interchange, which is preserved 
throughout the calculation. 

6 Conclusions 

Our explicit results are summarized in Table 1. The contributions to {g — 2)i of the 
graviton and gravitino sectors in differential renormalization are separately finite and 
independent of the scale M. In fact, in our procedure the logarithms of M keep track 
of the singularities and their sum is shown to vanish. Notice that these logarithms 
are proportional to the infinities of the regularization methods in momentum space, 
as expected at one loop. Obviously, finite parts may be different and are responsible 
for the preservation or the breaking of the symmetries in the problem. In differential 
renormalization the total contribution cancels out as dictated by supersymmetry. In 
the dimensional reduction scheme the total contribution also vanishes. This is not the 
case in dimensional regularization as this scheme is known to break supersymmetry. 

Note that the values of the {g — 2)i contribution of the graviton and gravitino 
sectors are different in differential renormalization (1,-1) and in dimensional reduc- 
tion (—1/2, 1/2), although both schemes preserve supersymmetry. These contributions 
appear to be regularization dependent. This ambiguity is related to the presence of 
linear singularities (divergences) at one loop in SUGRA. All that can be said is that 
global supersymmetry ensures a vanishing value of {g — 2)i and a cancellation of the 



ambiguities between the graviton and gravitino sectors. In Ref. [21| yet another finite 
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Diagram 



Differential Dimensional Dimensional 

Renormalization Reduction Regularization 



Dl 

D2+D3 

D4+D5 



ilog 



11 



4fiog 



M2 



25 1^^ 29 

18 3 ra-4 18 



n 11^^ 35 

18 3 n-4 9 



MS) 



+ 1 



-4;i^ + 6 



1^^ 61 

3 ra-4 36 



11 1 
3 n-4 



32 
9 



^n-4 ^ ' 



Graviton 
(D1+D2+D3+D4+D5) -1 



1/2 



7/4 



D6 

D7+D8 

D9+D10 



-llog(f^) 



+ 



19 8^^ 37 

18 3 n-4 18 



21og(|^) 



n-4 



+ 2 



8^^ 55 

3 n-4 18 



2 Ino- /'^^^ 4- 11 ^-k 1 4^^ 13 

3 ^'-^S V m^y "•" 18 3 n-4 9 3 n-4 9 



n-4 



+ 4 



Gravitino 
(D6+D7+D8+D9+D10) 1 



-1/2 



-1/2 



TOTAL 

(Graviton+Gravitino) 



5/4 



Table 1: Contributions of the diagrams in Fig. 2 to ( ^^ 1 in units of ^"^ , obtained 
with DR, dimensional reduction and dimensional regularization. 
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value for the graviton contribution to {g — 2)i, —13/12(Gm?/TT), was obtained using 
source theory techniques. The situation is that in four different schemes the graviton 
contribution is different. The Pauh term is ahowed in the gravity lagrangian and its 
one-loop corrections are scheme dependent in such a non-renormalizable theory. 

An important point developed in our work is the exact maintenance of symmetries 
within the differential renormalization scheme. Direct renormalization with different 
scales requires enforcing the relevant symmetries (gauge and supersymmetry) with the 
corresponding Ward identities at the end of the calculation. In our approach all ex- 
pressions have been handled in a symmetric way. We have reduced all singularities to 
a minimal set of independent functions. This was done relating diagrams of different 
topology with the technique of point separation. This procedure ensures that singu- 
larities are consistently treated. Then we obtain a result which is supersymmetric and 
compatible with U(l) gauge invariance. 
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Appendix A 



In this Appendix we give the Feynman rules (Figs, y and Q) for the lagrangian (| 



2.4). They are written in Euchdean space with the convention {70,7/3} = 25q,^. 

'^^ Yn gap 




-^np. 




-iCKda^ilp 



y^ 



ie{di - dl 






6l,r Sl,r 



^VL,R{f + ^hc, 



V2eVR,i 



eL,R 




■f7«(9^-9; 




fS ^f3 




i[^^7m]7o 



en 
/2 



'a'PL,R 







Figure 3: Feynman rules for vertices. Vr^l = 2(^ -'- 7^) ^'^^ ^^^ chiral projectors. All 
derivatives are with respect to the vertex space-time point. The superscripts indicate 
the field they are acting on. The rules for diagrams with opposite charge and fermion 
number arrows are obtained from these by the transformation FR -^ 75FRT75 [25]. 
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m)Ara{x-y) 



Am(x - y) 



y 



H / V 



y vAAAAAAA/^ X "/^'^ 



5^uMx - y) 



ivw^ X - ^''A(x - y) 



y vAAAAAAA/' X 



^IV 



LI V & p O 

AAA/v\A/>x 



2 [y^ipOva \ UpujOyp 

- 5^u5p„)l\{x - y) 



y ,^A,A,A,gA,A,A;A x l7- ^ Tm^Ix " 2/) 



Figure 4: Feynman rules for propagators. A(x) = 4;;;7-3- and Am(x) = __^ ™- u™^; g^j,g 
the massless and massive scalar propagators. Here the superscript on the derivatives 
refers to the space-time point. 



Appendix B 

In this Appendix we show how Ti[5Qn] and Ti[{dadp — ^J^^^nJD] can be expressed in 
terms of the logarithmically singular function Ti[n]. Let us first deal with the linearly 
singular Ti[(?an]: 

Ti[5«n] = A„(x)A„(y)aSn^-A(x-y) 

= -Amix)Amiy)d^Six-y) 

= -aS[A„(x)A„(y)<5(x - y)] + [(a>„(x))A„,(y)]<5(x - y) 

= 9^Ti[n] + P^A„,(x))A„(x)]<5(x-y). (B.l) 
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On the other hand, 



= -a^Ti[n]-[A„(x)(aSA„(x))]5(x-y). (B.2) 

So, combining both expressions 

Ti[dM = \{dl-dl)T,[U]. (B.3) 

With the 'integration by parts' prescription the hnear singularity has been reduced to 
a logarithmic one. For Ti[{dadp — ^6apC\)U\ we have 

T,[{do.dp-\5^pU)n] = -1(9-9- + a^9|-l5„^(n- + n^))[(A„(x))25(x-y)] 

+ [5-A„(x)a^A^(x) - i5„^a-A^(x)5-A„(x)] 

X 6{x - y) , (B.4) 

where the previous result has been used. Using recurrence relations among modified 
Bessel functions, the second term in the r.h.s. can be written 



1 .XaX[3 1 m^K^{mx) 



(4^2)2 V 3,2 4"«P^ ^2 

+ rn^{KQ[mx) — Kl{mx))] 

m?Kl(mx) m?KQ{mx)Ki{mx) 
X [5(rE - y){ -2 + 

+ m^{Kl{mx) - Kl{mx)))], (B.5) 

where the identity 

[d^.fixMx -y) = {dl + dl)[f{x)5{x - y)] (B.6) 

has been used to obtain the last equality. Only the term with ^ — - in Eq. ( p. 51 ) is 

singular and equal to 

- \Wa + di){d], + di) - \5^p{d; + dy){d; + dy)]T,p] . (B.7) 

Thus we finally obtain 
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Appendix C 

In order to extract the contribution to the anomalous magnetic moment of a given 
graph from its full renormalized expression, external fields must be on-shell. External 
derivatives in IPX graphs act directly by parts on the external fields, so the Dirac 
or Maxwell equations can be used straightforwardly. When internal derivatives are 
present, however, the situation is more involved due to the noncommutative character 
of the derivatives. A simple procedure to follow is to perform a Fourier transform of 
the whole graph and deal with momenta, which are commuting objects. Notice that 
all Fourier transforms are ultraviolet convergent since the singular pieces have already 
been renormalized. On-shell conditions can then be readily imposed and the {g — 2); 
contribution identified. 

The Fourier transform of a distribution f{x,y) is 

f{p,p') = fd'xd'ye'P--e'P'-yf{x,y) , (C.9) 



where p and p' are the incoming momenta of the external leptons and q = p + p' is the 
outgoing momentum of the external photon. External derivatives yield 

and we are left with the Fourier transforms of T-functions. When computing them we 
distinguish the 'originally regular' ones from the 'renormalized' ones (which, of course, 
are also regular, but involve different kinds of functions). The originally regular T- 
functions are products of three propagators with some inner derivatives. Their Fourier 
transform is a convolution of momentum space propagators with the corresponding 
internal momenta. For Ti-functions: 

fi[o{d)]=i[o{d^ik)], (c.ii) 

where 



l[fik)]^J j^A^{p + k)A^{p'-k)Aik) , (C.12) 

A(fe) = l, A„.(fc) = ^^. (C.13) 
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Thus, e.g., 



tl[9a9/3 - \Sa^O] = -l[kakp - ^kHafs]. (C.14) 



The integrals I [/(A;)] appear in one-loop momentum space calculations and can be 
evaluated with standard techniques such as Passarino-Veltman reduction and Feynman 
parametrization. In our case, since we are interested in the static limit 

p2=p'2 = -m2 , g2^0 , (C.15) 

it is convenient to impose these conditions before computing the integrals, which then 



become much simpler [26]. 



Formally, the Fourier transform of originally regular T2-functions is identical: 

T2[0{d)]=J[0{d^ik)] , (C.16) 

where 

J[/(fc)] ^ / ^ A(p + k)A{p' - k)AUk) . (C.17) 

However the presence of two massless propagators makes the static limit of these in- 
tegrals infrared divergent^, so one must keep g^ 7^ 0. It is gratifying to note that all 
the infrared divergencies cancel out in the {g — 2)/ part of each diagram, rendering the 
limit g^ — > finite. 

Renormalized expressions of T-functions always contain a delta function, 5{x — y). 
Their Fourier transforms become integrals over one (four-dimensional) variable: 



/ 



(fx d% e*P-^e*P ■y[gix)6ix - y)] = / d^x e''^-''g{x) , (C.18) 



with q^ = P^i -\- p'n.- Hence we only need the Fourier transforms of certain functions of 
one variable containing logarithms or modified Bessel functions. Derivatives, if present, 
are trivially pulled out yielding momenta q. The transform of logM'^x'^ / x"^ was given 
m||; 



Note that Eq. (C.19) is divergent when q^ -^ 0. In the computation of {g — 2)i, 
infrared divergencies coming from originally regular pieces and from renormalized pieces 
cancel separately. To obtain the Fourier transforms of expressions with modified Bessel 
functions, the following recurrence relations in four dimensions prove useful: 

□ ^iM = :^lM_4vr^5(,), (C.20) 



^Ti[l] is also infrared divergent when q^ — * 0, but it does not appear in the {g — 2); parts. 



2n 



aK2{x) = Ko{x)-8tt^5{x), (C.21) 

aKoix) = Ko{x)-2^^^, (C.22) 

X 

aKf{x) = 2{K^{x)+K'l{x))-ATr^6{x), (C.23) 

aiKfix) - Kiix)) = 4 ^o(^)-^i(^) _4^2^(^)_ (C.24) 

X 

All these expressions are tempered distributions. The transforms we need are 

^^^^,g.^mKi^rnx)_ ^ 47r^Am{q), (C.25) 

X 

f d^xe^'^-'^Koimx) = 8TT'^[Am{q)f, (C.26) 



™2 ^ j4;^ 



I d^x e^'^-^m^K'.imx) - K^{mx)) = -1287r4^ J —^ A„(g - fc)[A„(g)] 



47r2 
,2-.o 2 7r2 



+ .2 



- „ 2- (C-28) 

Finally, the momentum space expressions allow for direct use of the Dirac equation 
(^ -^ —im, 'p' -^ im), and the {g — 2)i parts are easily recognized, for they are pro- 
portional to Pfj_ — p'^. 
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